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ABSTRACT 

We  consider  the  problem  of  reconstructing  a  signal  field 
measured  by  a  large  scale  sensor  network  with  mobile  agents. 
Sensors  transmit  packets  containing  measurement  data  to 
mobile  agents  using  either  random  or  deterministic  medium 
access  control  (MAC)  schemes,  and  the  signal  field  is  recon¬ 
structed  by  mobile  agents  that  minimize  the  mean  square  er¬ 
ror  of  the  reconstruction.  For  the  one-dimensional  Gauss- 
Markov  field,  we  investigate  the  relation  between  the  sys¬ 
tem  throughput  and  reconstruction  distortion,  for  different 
types  of  MAC  schemes.  We  show  that  at  low  throughput 
level,  increasing  system  throughput  decreases  the  recon¬ 
struction  distortion  considerably.  But  the  improvement  is 
much  less  when  the  throughput  is  relatively  high.  We  also 
show  that  the  choice  of  MAC  schemes  can  affect  the  re¬ 
construction  performance  significantly,  especially  when  the 
measurement  noise  is  low. 

1.  INTRODUCTION 

A  reachback  sensor  network  with  mobile  agents  is  proposed 
in  [1]  as  a  new  network  structure  for  low  power,  large  scale 
sensor  networks.  In  such  network,  sensors  and  mobile  agents 
arc  two  types  of  operating  nodes.  Sensors  arc  densely  de¬ 
ployed  in  the  field  to  sense  information  or  take  measure¬ 
ments,  and  mobile  agents  arc  in  charge  of  collecting  infor¬ 
mation  from  sensors.  At  the  time  of  data  collection,  sensors 
transmit  their  data  directly  back  to  the  mobile  agents.  By 
fully  utilizing  mobile  agents’  capability,  such  sensor  net¬ 
work  structure  simplifies  node  level  processing  and  can  ben¬ 
efit  from  node  redundancy  for  energy  efficient  and  reliable 
communications. 

This  work  was  supported  in  part  by  the  Army  Research  Office  under 
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cation  and  Networks  under  Grant  DAAD 19-01 -2-001 1. 


One  type  of  applications  in  such  sensor  network  is  to 
provide  measurements  of  some  phenomena  of  interest  in  a 
field,  such  as  temperature.  Specifically,  the  physical  source 
information  in  a  field  A  can  be  viewed  as  a  continuous 
stochastic  process.  A  very  large  amount  of  sensors  are  spread 
in  the  field  for  local  measurements.  At  a  prearranged  time, 
each  sensor  senses  the  information  and  stores  its  measure¬ 
ment.  When  the  mobile  agent  is  ready  to  collect  data,  sen¬ 
sors  transmit  their  packets  back  to  the  mobile  agent  ac¬ 
cording  to  a  specified  medium  access  control  (MAC)  proto¬ 
col.  After  a  fixed  collecting  time,  the  mobile  agent  receives 
some  data  packets  from  some  locations  in  the  field.  Based 
on  these  samples,  we  reconstruct  the  original  signal  field. 

For  sensor  networks  with  random  access  MAC  schemes, 
due  to  the  access  randomness,  both  the  number  of  data  pack¬ 
ets  successfully  received  and  their  origins  arc  random.  This 
appears  as  if  the  mobile  agent  randomly  samples  (both  the 
location  and  the  number  of  samples)  the  measurement  data 
in  A  during  the  fixed  time  duration  T.  Intuitively,  the  more 
data  samples  collected,  the  better  signal  reconstruction  per¬ 
formance.  Notice  that  the  average  number  of  data  samples 
the  mobile  agent  receives  is  determined  by  the  through¬ 
put  of  the  network.  Then  how  does  the  reconstruction  per¬ 
formance  vary  with  the  throughput?  The  significance  of 
characterizing  the  relation  between  throughput  and  recon¬ 
struction  distortion  lies  in  two  folds:  first,  a  given  system 
throughput  provides  the  expected  field  reconstruction  per¬ 
formance  under  such  sensor  network;  second,  and  probably 
more  important,  to  target  the  reconstruction  at  a  specified 
distortion  tolerance  level,  it  indicates  how  much  through¬ 
put  the  network  is  required,  and  thus  directs  the  quality-of- 
service  (QoS)  specific  network  design.  Furthermore,  it  is 
also  important  to  know  the  variation  rate  of  distortion  un¬ 
der  different  throughput  levels,  as  it  indicates  the  effective¬ 
ness  of  throughput  improvement  on  the  reconstruction  per¬ 
formance.  Unlike  conventional  communication  networks 
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where  independent  data  arc  generated  from  each  user  nodes, 
having  (highly)  correlated  data  among  sensor  nodes  is  an 
important  feature  in  a  densely  deployed  sensor  network. 
Therefore,  another  related  question  to  ask  is  that  how  dif¬ 
ferent  types  of  MAC  schemes  affect  reconstruction  perfor¬ 
mance?  The  validity  of  such  question  is  evident  from  the 
fact  that  the  pattern  of  received  data  sample  locations  affects 
signal  reconstruction  performance.  Intuitively,  a  uniformly 
sampled  data  results  in  less  reconstruction  distortion  than 
that  from  using  data  samples  grouped  in  a  small  subarea. 

In  this  paper,  we  consider  the  problem  of  measuring  a 
one-dimensional  signal  field  using  a  reachback  sensor  net¬ 
work  with  a  mobile  agent.  The  source  signal  is  modeled  as 
a  Gauss-Markov  stochastic  process.  Based  on  the  received 
data  at  the  mobile  agent,  the  minimum  mean  square  error 
(MMSE)  estimator  is  used  for  signal  reconstruction.  The 
reconstruction  distortion  is  measured  by  the  maximum  esti¬ 
mation  error  of  reconstruction.  We  investigate  the  relation 
between  the  throughput  of  the  network  and  reconstruction 
distortion.  We  consider  both  noiseless  and  noisy  measure¬ 
ment  cases.  We  show  that  increasing  throughput  decreases 
the  distortion  considerably  at  low  throughput  level,  but  this 
effect  is  much  less  when  the  throughput  is  relatively  high. 
Besides  random  access  MAC  schemes,  we  also  evaluate  the 
performance  under  the  MAC  with  uniform  packet  recep¬ 
tion,  i.e.,  the  packets  received  are  from  uniformly  spaced 
locations  in  the  field.  Comparing  the  performance  under 
the  above  two  types  of  MAC  schemes,  we  show  that  the 
latter  type  of  MAC  schemes  significantly  improves  recon¬ 
struction  performance,  especially  at  low  measurement  noise 
level.  This  indicates  the  importance  of  sensor-network  spe¬ 
cific  MAC  designs  on  the  ultimate  system  performance. 

The  problems  on  sensor  network  communications  have 
attracted  a  growing  research  interest.  In  terms  of  medium 
access  control,  many  MAC  protocols  have  been  proposed 
aiming  to  the  special  needs  and  requirements  for  both  ad 
hoc  sensor  networks  [2-5]  and  reachback  sensor  networks 
with  mobile  agents  [6,7].  The  problem  of  signal  estima¬ 
tion  and  reconstruction,  admitting  the  specific  characteris¬ 
tics  under  such  sensor  network  settings,  is  still  a  relatively 
new  topic  for  investigation.  Perhaps  the  most  relevant  work 
is  [8],  where  the  problem  of  random  sampling  of  a  dynamic 
system  is  discussed,  and  sequential  estimation  of  a  dynamic 
source  is  considered  using  random  sampling  in  time. 

2.  PROBLEM  STATEMENT 
2.1.  Source  Model 

Consider  a  one-dimensional  field  of  length  D ,  denoted  by 
A  =  [0,  D\.  Let  SfX'x)  (x  €  A)  be  the  source  of  interest  in  A 


at  time  t.  We  assume  that  the  spatial  dynamic  of  Si(x)  is  a 
stationary  Gaussian  random  process  given  by  the  following 
lineal-  stochastic  differential  equation: 

dSt(x )  =  fSt{x)dx  +  adWt(x)  (1) 

where  /,  cr  €  R  are  known  and  /  <  01,  {Wt{x)  :  x  >  0} 

2 

is  a  standard  Brownian  motion,  and  St{x)  ~  AA(0,  jfj|)  . 

2.2.  Sensor  Measurement  and  Data  Reception  via  Mo¬ 
bile  Agent 

Suppose  that  a  very  large  amount  of  sensors  are  uniformly 
deployed  in  A,  and  every  sensor  knows  its  relative  location. 
The  sensor  measurement  at  time  t  is  given  by 

Yt(x)  =  St(x)+Nt(x)  (2) 

where  Nt(x)  is  an  independent  and  identically  distributed 
(i.i.d.)  zero  mean  white  Gaussian  measurement  noise  with 
variance  o^  .  Each  sensor  stores  its  local  measurement  of 
signal  and  waits  for  the  mobile  agent  to  collect.  Note  that 
we  will  reconstruct  the  source  signal  {St(x)  :  x  €  A}  at 
time  t  based  only  on  the  sensor  measurements  at  that  time: 
\Yf  {x)  :  x  €  A}.  Therefore,  we  drop  the  time  index  for 
brevity  in  the  following  presentation. 

When  the  mobile  agent  is  ready  for  data  collection,  sen¬ 
sors  transmit  packets,  consisting  measurement  data  along 
with  its  location  information,  to  the  mobile  agent  through 
a  common  wireless  channel,  using  either  random  or  deter¬ 
ministic  MAC  schemes.  To  remove  the  boundary  effect  for 
signal  reconstruction,  we  assume  that  the  mobile  agent  al¬ 
ways  first  obtains  packets  from  the  two  boundaries  of  A, 
i.e.,  T' (0)  and  Y  (D),  before  starting  the  collection  of  the 
rest  packets.  After  some  fixed  collecting  time  T,  assume 
that  the  mobile  agent  receives  total  M  packets  originated 
from  some  M  points  in  A.  Based  on  these  data  samples, 
we  reconstruct  the  signal  field  by  smoothing  technique  us¬ 
ing  the  MMSE  estimator.  Fig.  1  shows  an  example  of  a 
resulting  sampling  of  the  signal  field  in  A.  Given  the  re¬ 
ceived  data  at  points  p«  =  [0,  Pi,  •  •  •  ,  Pm,  D],  the  signal 
reconstruction  by  MMSE  smoothing  is  then  given  by 

S{x)  =  E{S(x)\Yp,p  e  pm},  x  €  A.  (3) 

If  a  random  access  MAC  scheme  is  used  in  the  sensor 
network,  both  the  number  of  data  packets  successfully  re¬ 
ceived  (M)  and  their  origins  (pM)  are  random.  This  ap¬ 
peal's  as  if  the  mobile  agent  randomly  samples  the  sensor 
measurement  data  in  A  (both  the  location  and  the  number 

1  f  <  0,  thus,  (1)  admits  a  stationary  solution. 
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of  samples)  during  collection  time  T.  Let  A  be  the  through¬ 
put  of  the  network.  We  assume  that  the  number  of  received 
packets  M  during  time  T  is  Poisson  distributed  with  rate  A: 


P(M  =  n ) 


=  e 


-AT 


(A  T)n 
n\ 


Then,  the  average  number  of  received  data  samples  during 
T  is  given  by  E{M}  =  XT. 


S(x) 


D 


Fig.  1:  a  1-D  signal  field  sampled  by  a  sensor  network 


2.3.  Source  Reconstruction  Distortion 

Given  M  +  2  received  packets  from  pM,  we  define  the  field 
reconstruction  distortion  by  the  maximum  estimation  error 
in  A: 

£(Pm)  =  maJE{\S(x)  -  S{x)\2  \pM}.  (4) 

Under  throughput  A  and  collecting  time  T,  the  average  dis¬ 
tortion  is  then  given  by 

£(XT)  i  E{£(  pM)}  (5) 

=  E{E{£( pM)  |  M}}  (6) 

where  in  (6),  the  second  expectation  is  taken  over  pM  for  a 
given  M,  and  the  first  one  is  taken  over  M. 

Clearly,  the  average  distortion  £{XT)  is  a  function  of 
the  throughput  A  of  this  sensor  network.  Intuitively,  the  net¬ 
work  with  higher  throughput  results  in,  on  average,  more 
data  samples  received  at  the  mobile  agent,  which  provides 
better  estimate  of  the  signal  field  and  results  in  lower  recon¬ 
struction  distortion.  In  the  following,  we  analyze  how  the 
reconstruction  distortion  5  (AT)  varies  with  the  throughput 
A. 


3.  SOURCE  RECONSTRUCTION:  DISTORTION 
VS.  THROUGHPUT 

We  assume  that  under  the  random  access  MAC  scheme,  the 
data  packet  generated  from  any  sensor  has  equal  chance  to 
be  successively  received  by  the  mobile  agent.  Therefore, 
given  there  arc  M  packets  received2 ,  the  packet  generation 

2  The  actual  number  of  received  packets  is  M  +  2.  For  convenience, 
we  only  count  those  packets  not  from  the  two  boundaries  of  A. 


locations  {Pi,---,  Pm}  are  i.i.d.  random  variables  with 
uniform  distribution  U( 0,  D). 

To  contrast  with  such  random  access  MAC  schemes,  we 
also  consider  another  type  of  MAC  schemes,  referred  to  as 
the  deterministic  MAC  with  uniform  reception.  This  type 
of  MAC  schemes  is  designed  in  such  a  way  that,  instead 
of  their  generation  locations  being  random,  the  M  received 
packets  at  the  mobile  agent  arc  originated  from  uniformly 
spaced  locations  in  A,  i.e.,  pM  =  {0,  ,  ■  ■  ■  >D}. 

In  the  following,  we  evaluate  the  average  distortion  un¬ 
der  both  types  of  MAC  schemes,  and  compare  their  perfor¬ 
mance. 


3.1.  Sensors  with  Noiseless  Measurement  Data 


We  first  consider  the  case  when  sensor  measurements  con¬ 
tain  no  noise,  i.e.,  Y(x)  =  S(x).  Given  M  =  m  received 
packets,  the  corresponding  order  statistics  of  the  packet  gen¬ 
eration  locations  {Pi,  •  ■  ■ ,  Pm},  denoted  by  P( i)  <  •••  < 
P(rn) ,  have  the  joint  density  given  by 

r  /  1  _  /  15^  if  0  <  Xl  <  •••  <  Xm  <  D 

JPw,-,P(m)\x u  ,xm)  -  otherwise. 

(V) 

From  (1),  due  to  the  Markov  property  of  the  stochastic  pro¬ 
cess,  given  samples  at  locations  andp^+1^,  the  signal 
S(x),  where  plt  >  <  x  <  Py+y,  is  completed  determined 
by  S(p^y)  and  S(p^i+y).  By  the  Gaussian  property  of  the 
process  {5(a:)},  the  MMSE  estimator  in  (3)  is  then  given 
by 

S(x)  =  E{S(x)S^H}E-1{S^S^H}S^  (8) 

forp(.)  <  X  <  p{i+1y  Where  =  [5(p(q),  S(p(i+i))]T. 

Let 

di  =  P(j+i)  -  P(i)  (9) 

for  i  =  0,  •  •  •  ,  M?  Then,  the  corresponding  maximum 
MMSE  of  S(x),  for  p^y  <  x  <  Py+y,  is  obtained  at  the 
middle  point  of  pyy  and  P(*+i),  and  it  is  given  by 


_<S(a:)l2l  P(0’P(U1)} 

=  max  {P|S'(a:)|2  — 

P(i)<*<P(i  + 1) 


S{5(*)S«H}S-1{S«sW"}£;ff{5(a:)S«ff}} 

1  —  efdi  a2 

1  +  efdi  1 2/ 1  ^ 


Therefore,  the  distortion  £  (pM)  in  (4)  is  determined  by  the 
maximum  of  distances  between  any  two  consecutive  data 

3we  denote  P(0)  =  0  and  P(m+i)  =  D. 
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samples 


For  deterministic  MAC  with  uniform  reception,  £(XT)  can 
be  readily  obtained  as  follow 


£{dmax)  —  £(Pm) 


^  max  jj2 

1  +  e^'c!max  | 2/ | 


(11) 


where 


£eq(AT)  =  e 


-AT 


<7 

M 


E 


f-tL. 

g*7  m+1 


m=0  V 1  +  e  ™+! 


(AT)r 


m\ 


(14) 


((max  —  max  (Lj . 

0<i<M 

Thus,  to  calculate  the  average  distortion  £(XT),  we  only 
need  to  find  the  probability  distribution  of  the  maximum 
sample  distance. 

Under  the  random  access  MAC  schemes,  given  M  = 
m,  we  now  find  the  cdf  of  the  maximum  sample  distance 
F<liw<y  (-|m).  By  the  transformation  of  random  variables, 
from  (7)  and  (9),  we  have  the  joint  distribution  of  sample 
distances  {do,  ■  ■  •  ,  dm~ i}  given  by 


/do,  -  ,dm-i  (®0j  ->Xm—  l) 


{ml 

r 


if  Xi  >  0,  Yh= o1  Xi<D 
otherwise. 


and 


Fig.  2  shows  the  curves  of  throughput-distortion  char¬ 
acteristics  under  the  two  types  of  MAC  schemes.  We  set 
/  =  —0.1,  cr2  =  0.2,  and  T  =  100.  We  observe  that  when 
A  is  low,  the  average  distortion  decreases  fast  with  increas¬ 
ing  throughput  A.  However,  as  A  becomes  larger,  the  de¬ 
creasing  rate  of  distortion  becomes  slower.  In  other  word, 
the  signal  reconstruction  performance  is  more  sensitive  to 
the  throughput  improvement  at  low  throughput  level.  Com¬ 
paring  the  performance  under  random  MAC  schemes  with 
that  under  the  deterministic  MAC  scheme  with  uniform  re¬ 
ception,  we  observe  that  a  significant  performance  gain  can 
be  achieved  using  the  latter  scheme.  This  indicates  that, 
in  such  sensor  networks  with  highly  correlated  data  infor¬ 
mation,  the  MAC  design  with  uniform  reception  is  highly 
desired. 


m—1 


^dmax  ix\m)  =  Pr(do  <  X,  ■  ■  ■  ,  dm_l  <  x,  D  -  ^2  di  <  x) 


HX 

g{x  or 


i= 0 


i  xm—\)dxQ  •  •  •  dxm— i (12) 


where 


g{x  o,  ,  xm— i) 

—  /do,-"  ,dm-i  (^Oj  5  Xm-1)^[D~YX£o  Xi<x]  (^O)  *xm-l) 

and  1  b(x)  is  the  indicator  function 

1,  if  x  E  B 


1  b{x)  = 


0,  if  x  ^  B. 


T=100,  o2=1 ,  f=-0.1 ,  Noiseless 
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Notice  that  to  calculate  Pdmax  (xjm)  in  (12)  is  essentially  to 
calculate  the  volume  of  a  m-dimensional  hypercube  under 
the  two  hyperplane  constraints.  The  resulting  formula  for 
i'dmax  (x\m)  is  given  in  (17).  Thus,  the  average  distortion  of 
signal  reconstruction  £(XT)  in  (6)  is  obtained  by 


00  /  fOQ 

£(XT)  =  [  Jo  P(£(dmax)  >  x ) 


—  AT 


m= 0 
oo 

E 

m= 0  L 


dx  )  e 


-xt  (AT)T 


ml 


/  TO))  dx 

(AT)™ 


-\-£ 


D 
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Fig.  2:  Distortion  vs.  throughput 

3.2.  Sensors  with  Noisy  Measurement  Data 

When  noises  arc  presented  in  the  signal  measurements,  all 
received  data  samples  need  to  be  used  to  reconstruct  signal 
at  each  point  in  A.  Let  SPM  =  [5(P(0)),  •  •  •  ,  S,(P(M))]T 
be  the  source  signal  of  which  the  noisy  measurements  arc 
received  at  the  mobile  agent.  In  this  case,  the  MMSE  of 
S(x)  is  given  by 

M{x-,pM)=E{\S{x)  -  S,(x)|2  |  pM} 
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where  r(x'.  pM)  =  ^{5(x)S^m}  and  RPm  =  E{SPmS^m}. 

Hence,  £  (AT)  is  obtained  as 

°°  (  fD  rD 

£{\T)=e~XT  ^  /  ■  ■  /  maxM(j/;xm) 

TO=0  V°  ^  " 

(XT)171 

,P(m)  (x1j  ' ' ' ,*m)  dx1---dxm)  V  —  .  (15) 

The  distortion  under  the  MAC  with  uniform  reception  can 
be  obtained  similarly 

£eq(AT)  =  e“AT  M(z;  pj^-  (16) 

m=0 

where  in  this  case  pm  =  [0,  •  ,D]. 

Fig.  3  shows  the  distortion  curve  under  the  variation  of 
throughput  at  SNR  =  HMB  and  20dB,  respectively.  The 
same  setup  as  in  Fig.  2  is  used.  We  observe  that  as  the  noise 
becomes  higher,  the  curve  become  more  flat.  This  indicates 
that  due  to  the  measurement  noise,  throughput  improvement 
is  less  effective  on  the  reconstruction  performance  gain.  Fur¬ 
thermore,  we  observe  that  the  performance  gap  under  the 
two  different  MAC  schemes  becomes  smaller  as  the  mea¬ 
surement  noise  become  higher. 


T=100,  o2=1,  f=-0.1,  SNR=20dB 


Fig.  3  :  Distortion  vs.  throughput 


4.  CONCLUSION 

In  this  paper,  the  problem  of  reconstructing  a  signal  field 
measured  by  a  densely  deployed  sensor  network  with  a  mo¬ 
bile  agent  is  considered.  We  investigate  the  relation  be¬ 
tween  source  reconstruction  performance  and  network  through¬ 
put  under  both  random  and  deterministic  MAC  designs.  We 
show  that  throughput  improvement  is  more  effective  on  re¬ 
construction  performance  gain  at  a  low  throughput  level 
than  that  at  a  relatively  high  level.  Furthermore,  we  show 
that  designing  a  MAC  scheme  with  uniform  reception  can 
obtain  a  considerable  gain  on  reconstruction  performance. 

5.  REFERENCES 

[1]  L.  Tong,  Q.  Zhao,  and  S.  Adireddy,  “Sensor  networks  with 
mobile  agents,”  in  Proc.  IEEE  MILCOM’03 ,  Boston,  MA, 
October  2003. 

[2]  R.  Iyer  and  L.  Kleinrock,  “QoS  Control  for  Sensor  Net¬ 
works,”  in  Proc.  ICC  2003,  May  2003. 

[3]  K.  Sohrabi,  J.  Gao,  V.  Ailawadhi,  and  G.  Pottie,  “Protocols 
for  self-organization  of  a  wireless  sensor  network,”  IEEE 
Personal  Communication,  pages  16-27,  October  2000. 

[4]  A.  Woo  and  D.  Culler,  “A  transmission  control  scheme  for 
medium  access  in  sensor  networks,”  in  Proc.  ACM/IEEE 
Conference  on  Mobile  Communication  and  Networking, 
pages  221-235,  Rome,  Italy,  July  2001. 

[5]  W.  Ye,  J.  Heidemann,  and  D.  Estrin,  “An  energy-efficient 
MAC  protocols  for  wireless  sensor  networks,”  in  Proc.  of 
IEEE  INFOCOM,  pages  1567-1576,  New  York,  June  2002. 

[6]  S.  Adireddy,  P.  Venkitasubramaniam,  and  L.  Tong,  “Oppor¬ 
tunistic  ALOHA  and  cross  layer  design  for  sensor  networks,” 
in  Proc.  IEEE  MILCOM’03,  Boston,  MA,  October  2003. 

[7]  Q.  Zhao  and  L.  Tong,  “QoS  specific  information  retrieval 
for  densely  deployed  sensor  network,”  in  Proc.  IEEE  MIL- 
COM’03,  October  2003. 

[8]  M.  Micheli  and  M.  I.  Jordan,  “Random  Sampling  of  a 
Continuous-time  Stochastic  Dynamical  System,”  Proc.  of  the 
15th  International  Symposium  on  the  Mathematical  Theory 
of  Networks  and  Systems  (MTNS  2002),  University  of  Notre 
Dame,  South  Bend,  Indiana,  August  2002. 


5 


